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Abstract 

The first observation of (anti)deuterons in deep inelastic scattering at HERA lias 
been made with the ZEUS detector at a centre-of-mass energy of 300-318 GeV 
using an integrated luminosity of 120 pb~^. The measurement was performed 
in the central rapidity region for transverse momentum per unit of mass in the 
range 0.3 < px/M < 0.7. The particle rates have been extracted and inter- 
preted in terms of the coalescence model. The (anti)deuteron production yield is 
smaller than the (anti)proton yield by approximately three orders of magnitude, 
consistent with the world measurements. 
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1 Introduction 



Light stable nuclei, such as deuterons (d) and tritons (t), are loosely bound states whose 
production mechanism in high-energy collisions is poorly understood. Most measurements 
of light stable nuclei have been performed for antideuterons {d). A selection of d from 
primary interactions is more difficult as it requires separation of such states from particles 
produced by interactions of colliding beams with residual gas in the beam pipe and by 
secondary interactions in detector material. The first observation of d [1] was followed 
by a number of experiments on antideuteron production. The production rate of d in 
e+e~ — > qq collisions [2-5] is significantly lower than that measured in T(IS') and T(25') 
decays [2,5]. The d rate in e~^e~ —>■ qq is also lower than that in proton-nucleus (pA) [6,7], 
proton-proton (pp) [8] and photon-proton (7^) collisions at HERA [9], but higher than that 
in nucleus-nucleus collisions [10,11]. For heavy-ion collisions, the coalescence model [12] 
was proposed to explain the production of d{d). 

This paper presents the results of the first measurement of d and d in the central rapidity 
region of deep inelastic ep scattering (DIS). The analysis was performed for exchanged 
photon virtuality, Q"^, above 1 GeV^. 



2 Coalescence model for (anti)deuteron formation 

According to the coalescence model [12] developed for heavy-ion collisions, the production 
rate of d is determined by the overlap between the wave-function of a proton (p) and a 
neutron (n) with the wave-function of a d. In this case, the d cross section is the product of 
single-particle cross sections for protons and neutrons, with a coefficient of proportionality 
refiecting the spatial size of the fragmentation region emitting the particles. The same 
approach applies for d production. This model was also used to describe d{d) production 
in pp [8], 7p [9] and e"'"e~ [2,4] interactions. 

Assuming that all baryons are uncorrelated and the invariant differential cross section for 
neutrons is equal to that for protons, the invariant differential cross section for deuteron 
formation can be parameterised as 



where Ed{p) and ad{p) are the energy and the production cross section of the d(p), respec- 
tively, Pd{Pp) is the momentum of the d{p) and atot is the total ep cross section for the 
considered kinematic range. The coalescence parameter, B2, is inversely proportional to 
the volume of the fragmentation region emitting the particles. The same relation holds 
for d and p. If B2 is the same for particles and antiparticles, then the production ratio 
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d/d is equal to {p/pY- The coalescence parameter can be obtained from 




where M^^p) is the mass of the d{p), 7^ = Ea/Ma, R{d/p) is the ratio of the number of d 
to p expressed as a function of px/Mai^p-^, with px being the transverse momentum [9]. 



3 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [13]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [14], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar-angl^ region 15° < ^ < 164°. The transverse-momentum resolution for full-length 
tracks is a{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pr, with pr in GeV. To estimate the 
ionisation energy loss per unit length, dE/dx, of particles in the CTD [15], the truncated 
mean of the anode-wire pulse heights was calculated, which removes the lowest 10% and 
at least the highest 30% depending on the number of saturated hits. The measured 
dE/dx values were corrected by normalising to the average dE/dx for tracks around the 
region of minimum ionisation for pions with momentum p satisfying 0.3 < p < 0.4 GeV. 
Henceforth, dE/dx is quoted in units of minimum ionising particles (mips). The resolution 
of the dE/dx measurement for full-length tracks is about 9%. 

The high-resolution uranium-scintillator calorimeter (CAL) [16] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/ \^ for electrons 
and a{E)/E = 0.35/V^ for hadrons, with E in GeV. A presampler [17] mounted in front 
of the calorimeter and a scintillator-strip detector (SRTD) [18] were used to correct the 
energy of the scattered electroij^. The position of electrons scattered close to the electron 
beam direction is determined by the SRTD detector. 

The inactive material between the interaction region and the CTD, relevant for this 
analysis, consists of the central beam pipe made of aluminum with 1.5 mm wall thickness 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 

^ Henceforth the term electron is used to refer both to electrons and positrons. 
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and the inner diameter of 135 mm. The CTD inner wall with a diameter of 324 mm 
consists of two aluminum skins, each 0.7 mm thick, separated by a 8.6 mm gap filled with 
polyurethane foam with a nominal density of 0.05 g/cm^. 

The luminosity was measured using the bremsstrahlung process ep — > epy with the 
luminosity monitor [19], a lead-scintillator calorimeter placed in the HERA tunnel at 
Z = -107 m. 

4 Monte Carlo simulation 

To study the detector response, the Ariadne 4.12 Monte Carlo (MC) model [20] for 
the description of inclusive DIS events was used. The Ariadne program uses the Lund 
string model [21] for hadronisation, as implemented in Pythia 6.2 [22-24]. In its original 
version, this MC does not include a mechanism for the production of d or other light 
stable nuclei. To determine reconstruction efficiencies, a second Ariadne sample was 
generated in which (i's were included at the generator level by combining p and n with 
similar momenta. 

The Ariadne events were passed through a full simulation of the detector using the 
Geant 3.13 [25] program. The Geant simulation uses the Gheisha model [26] to 
simulate hadronic interactions in the material. The Geant program cannot be used for 
d as this particle is not included in the particle table. 

5 Event sample 

5.1 DIS event selection 

The data sample corresponds to an integrated luminosity of 120.3 pb"^ taken between 
1996 and 2000 with the ZEUS detector at HERA. This sample consists of 38.6 pb~^ of 
e~^p data taken at a centre-of-mass energy of 300 GeV, 65.0 pb~^ taken at 318 GeV and 
16.7 pb~^ of e~p data taken at 318 GeV. 

The search was performed using DIS events with exchanged-photon virtuality > 
1 GeV^. The event selection was similar to that used in a previous ZEUS publication [27]. 
A three-level trigger [13] was used to select events online. At the third-level trigger, an 
electron with an energy greater than 4 GeV was required. Data below ^ 20 GeV^ 
were prescaled to reduce trigger rates. 

The Bjorken scaling variable, xbj, and were reconstructed using the electron method 
(denoted by the subscript e), which uses measurements of the energy and angle of the 
scattered electron. The scattered-electron candidate was identified from the pattern of 
energy deposits in the CAL [28] . In addition, the inelasticity was reconstructed using the 
Jacquet-Blondel method [29], yjB, or the electron method, Ue- 
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For the final DIS sample, the following requirements were imposed: 

• Ql>l GeV^; 

• the impact point of the scattered electron on the RCAL outside the (X, Y) region 
(±12, ±6) cm centred on the beamline; 

• E^i > 8.5 GeV, where E^i is the energy of the scattered electron measured in the CAL 
and corrected for energy losses; 

• 35 < 5 < 65 GeV, where 5 = X] -^'(1 "^os 6'j), Ei is the energy of the i-th calorimeter 
cell, 9i is its polar angle and the sum runs over all cells; 

• Ve < 0.95 and yjB > 0.01; 

• at least three tracks fitted to the primary vertex to ensure a good reconstruction of 
the primary vertex and to reduce contributions from non-ep events; 

• I Zvtx |< 40 cm and \J X'^^^ + Y^^^ < 1 cm, where Z^tx, -^vtx and Yvtx are the coordinates 
of the vertex position determined from the tracks. 

The average of the selected sample was about 10 GeV^. 

5.2 Track selection and the dE/dx measurement 

The present analysis is based on charged tracks measured in the CTD. The tracks were 
required to have: 

• at least 40 CTD hits, with at least 8 of them for the dE/dx measurement; 

• the transverse momentum > 0.15 GeV. 

These cuts selected a region where the CTD track acceptance, as well as the resolutions 
in momentum and the dE/dx, were high. 

To identify particles originating from ep collisions, the following additional variables were 
reconstructed for each track: 

• the distance, AZ, of the Z-component of the track helix to Z^tx] 

• the distance of closest approach {DC A) of the track to the beam-spot location in the 
transverse plane. The beam-spot position is determined from the average primary- 
vertex distributions in X and Y for each data-taking period. The DC A is assigned a 
positive (negative) value if the beam spot lies left (right) of the particle path. 

Figure [1] shows the dE/dx distribution as a function of the track momentum for positive 
and negative tracks. The events were selected by requiring at least one track with dE/dx > 
2.5 mips. To reduce the fraction of tracks coming from non-ep collisions, the tracks were 
required to have |AZ| < 1cm and |DCy4| < 0.5 cm. After such a selection, clear bands 
corresponding to charged kaons, protons and deuterons were observed. The requirement 
dE/dx > 2.5 mips enhances the fraction of events with at least one particle with a mass 
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larger than the pion mass and leads to the discontinuity near dE/dx = 2.5 mips seen in 
Figure [H The lines show the most probable energy loss calculated from the Bethe-Bloch 
formula [30]. The dE/dx bands for K~ and p are slightly shifted with respect to the 
Bethe-Bloch expectations due to the geometrical structure of the CTD drift cells which 
leads to a different response to negative and positive tracks. 

Figure [2] shows the reconstructed masses, M, for different particle species. The masses 
were calculated from the measured track momentum and energy loss using the Bethe- 
Bloch formula. The mass distributions were fitted with asymmetrical Gaussian functions. 
The relative width obtained was 11% (7%) for the left (right) part of the function. 

The number of p{p) candidates in the mass region 0.7(0.6) < M < 1.5 GeV was 1.61 x 10^ 
(1.66 X 10^). Due to a shift in the dE/dx for negative tracks, the lower mass cut for p 
was at 0.6 GeV. The numbers of d and d in the mass window 1.5 < M < 2.5 GeV were 
309 and 62, respectively. The number of p migrating to the d mass region was estimated 
to be less than 1% of the total number of d candidates. A similar estimate was obtained 
for antiparticles. A small number of triton candidates was observed in the mass window 
2.5 < M < 3.5 GeV. However, due to low statistics, it was difficult to establish a peak 
inside this mass window, therefore, no conclusive statement on the origin of the tracks in 
the region 2.5 < M < 3.5 GeV was possible. 

The observed p{p) and d{d) candidates were required to be in the central rapidity region, 
\y\ < 0.4, and to have 0.3 < px/M < 0.7. This determines the kinematic range used for 
the cross-section calculations. 

5.3 Identification of particles produced in ep collisions 

The observed p{p) and d{d) candidates selected after the dE/dx mass cuts can originate 
from secondary interactions in the inactive material between the interaction point and 
the central tracking detector. 

In order to select p{p) and d{d) originating from ep collisions, both DC A and AZ cuts 
were removed and a statistical background subtraction based on the DC A distribution 
was performed. The AZ distributions for p{p) and d{d) after the mass cuts are shown in 
Figure [31 Clear peaks at AZ = are observed. To optimize the signal-over-background 
ratio for the DC A distribution, all candidates were selected using the | AZ \< 2(1) cm 
restriction for p,p [d, d). 

Figure H] shows the DC A distributions for p{p) and d{d) candidates. The distributions 
show peaks at zero due to tracks originating from the primary vertex. The number 
of particles originating from primary ep collisions was determined using the side-band 
background subtraction. A linear fit to the DC A distribution on either side of the peak 
region in the range 2 <| DC A \< 4 cm was performed. Then, the expected number of 

^ An asymmetric Gaussian has different widths for the left and right parts of the function. 
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background events in the signal region of | DC A \< 1.5(0.5) cm for p,p {d,d) candidates 
was subtracted. 

The number of p{p) obtained after the DC A side-band background subtraction was 1.52 x 
10^ (1.62 X 10^) . The numbers of d and d particles were 177 ± 17 and 53 ± 7, respectively. 
The difference in the observed numbers of p and p can be explained by different dE/dx 
efficiencies and the mass cuts for positive and negative tracks. Such a difference in the 
efficiencies for particles and antiparticles cannot explain the difference in the observed 
numbers of d and d. 

Figure [5] shows the distributions for several DIS kinematic variables: Xe, E^i and 5. In 
addition, rapidity {y) distributions for the selected candidates are shown. The numbers of 
p{p) and d{d) candidates were calculated in each bin from the DC A distributions after the 
side-band background subtraction. The distributions for d are consistent with those for p 
and p, while the d sample shows some deviations for the E^t variable and, consequently, 
for the 5 variable. 

6 Studies of background processes 

The following two background sources for heavy stable charged particles were considered: 

• interactions of the proton (or electron) beam with residual gas in the beam pipe, 
termed beam-gas interactions; 

• secondary interactions of particles in inactive material between the interaction point 
and the central tracking detector. 

6.1 Beam- gas interactions 

The contribution from proton-gas interactions is significantly reduced after the ZEUS 
three-level trigger which requires a scattered electron in the CAL. In addition, the re- 
quirement to accept only events with more than three tracks fitted to the primary vertex 
significantly diminishes the contribution from both electron-gas and proton-gas events. 
The remaining fraction of beam-gas interactions can be assessed by studying the Z^tx 
distribution. 

Figure [6] shows the Z^tx distributions for events with at least one p{p) or d{d) candidate. 
The distributions were reconstructed in the signal region | AZ |< 2(1) cm and | DC A |< 
1.5(0.5) cm for p,p {d,d) candidates without the background subtraction. Figure [6] shows 
that there is essentially no beam-gas background for d events. A small background for 
d at positive Z^t^ is expected from the DIS MC generated for inclusive DIS events in 
which (i's are solely produced by secondary interactions in the material in front of the 
CTD. This background is expected to have a fiat DC A and, therefore, is subtracted by 
the procedure described in Section 15. 3[ 
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The Zvtx distributions were fitted using a Gaussian function witli a first-order polynomial 
for the background description. The extracted Gaussian widths are fully consistent with 
those obtained for inclusive DIS events without the d preselection. 

To further study the Zvtx distribution, a special event selection was performed for non- 
colliding electron and proton bunches. Since the requirement to detect an electron with 
energy E^i > 8.5 GeV significantly reduces the rate of such background events, this re- 
quirement was not applied. All other tracking cuts were the same as in the d{d) selection. 
The requirement to accept events with at least three tracks fitted to the primary ver- 
tex rejects most of the beam-gas events (~ 95% from the total number of the triggered 
events). As expected, the remaining events show clear peaks at zero for the AZ and DC A 
distributions, but the reconstructed Z^tx distribution did not show a peak at zero. 

The enhancement at large Z^t^ for d, which was found to be consistent with that orig- 
inating from secondary interactions, could partially be due to electron-gas interactions. 
If one assumes that the background seen in Fig. MJo) is due to non-ep interactions, then 
the contribution from beam-gas interactions does not exceed 17% of the total number of 
events with a deuteron. 

6.2 Secondary interactions on inactive material 

A pure sample of DIS events will still contain deuterons produced by secondary interac- 
tions of particles in material. The aim of the side-band background subtraction discussed 
in Sect. I5.3l was to remove such a background contribution, assuming that the background 
processes do not create a residual peak at AZ = and DC A = 0. Several checks of this 
assumption are discussed below. 

The DC A and AZ distributions were investigated using a MC simulation of inclusive 
DIS events without d{d) production at the generator level. Deuterons from secondary 
interactions were selected as for the data. The reconstructed DC A and AZ for d did 
not show a peak at zero. A more detailed study of the DC A and AZ distributions was 
possible for p not originating from an ep collision at the MC generator level, since in 
this case the available MC statistics is significantly higher than for the d case. After the 
track-quality cuts, no peak at zero was observed in the DC A and AZ distributions. 

If a deuteron is produced by secondary interactions of the particles from the DIS event in 
the surrounding matter, the secondary d will not point precisely back to the interaction 
point, and both DC A and AZ distributions will be wider than in case of d and p. There- 
fore, the DCA and AZ distributions were fitted with double- Gaussian distributions to 
establish the width of the distributions. It was found that the observed deuteron DCA 
and AZ widths were consistent with the corresponding widths for p and p. 

One possible source for d is the reaction N + N — > d + ir, where one of the nucleons 
originates from ep collision, while the other one originates from the detector material in 
front of the CTD. For low initial nucleon momenta, the DCA of the d track is in general 
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large and it does not form an important background; at high initial nucleon momenta 
however, the DC A can become small enough that misidentification could become im- 
portanlo. Since the processes N + N ^ d + n can lead to an additional charged pion, 
this source of background deuterons can be studied by comparing the average charged 
multiplicity of tracks for d and d events. In addition, the distance of closest approach, 
DCA12, between the d track and other non-primary tracks in the same event should have 
an enhancement at zero. The study indicated that the average number of tracks for d 
events is smaller than that for d events. The rejection of events with |Z)CA12| < 2 cm did 
not lead to a statistically significant reduction in the number of the observed d events. 

Secondary deuterons may also be produced in pickup {p + n ^ d) reactions by primary 
p{n) interacting in the surrounding material. These deuterons, peaking in the direction 
of the primary p{n), point approximately to the interaction point and are therefore a 
potentially dangerous source of background. Experimental data on the pickup reactions 
at the relevant energy are scarce and therefore only a rough estimate of the size of this 
background is possible. From the extrapolation of data on Sm^^^ [31] and C [32, 33] 
targets using the K. Kikuchi theory [34] to allow for the change of material, the estimated 
d background from the pickup reaction was in the range 1 — 10% of the total number of 
observed d events, depending on the extrapolation input. 

The angular distributions of d from pickup reactions have also been investigated in several 
experiments [32,35,36] for various targets and for a range of p/M similar to the present 
analysis. In all cases, the angular distribution of d observed in these experiments would 
lead to a much wider DC A than that shown in Figure |l](b). 

7 Detector corrections 

In this analysis, all measurements are based on event ratios, therefore, the detector correc- 
tions due to DIS event selection and trigger efficiency were found to be small and thus are 
not discussed here. The detector corrections for the tracking efficiency and the efficiency 
of the dE/dx cuts are described below. 

7.1 Tracking efficiency 

The efficiency due to the track reconstruction, e, was estimated separately for p (p) and 
d using the Ariadne MC model (with d included at the generator level). The obtained 
efficiencies are about 0.95 for p and d and 0.90 for p. 

The method cannot be applied to d which are not treated in the Geant simulation. 
Therefore, the tracking efficiency for d was modelled as e{d) = e{d)e{p)/e{p). In the 
expression above, the hit reconstruction efficiency is described by the first term, e{d), 

^ Note that the cross section for the reaction N + N ^ d + tt decreases rapidly with increasing energy. 
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while the absorption loss (including annihilation) of d and p are assumed to be similar. 
This modelling assumes that the cross sections of annihilation in the detector material 
are the same for d and p, since the inelastic nuclear cross section of p is much larger than 
that of n for the momentum region less than 0.4 GeV [37]. The use of the geometrical 
model discussed in [11,37] and the model in which the p and n inelastic absorption cross 
sections are added linearly [4,37] to obtain the inelastic nuclear cross section of d, reduces 
e{d) by 1% and 5%, respectively. 

7.2 Efficiency of the dE/dx cuts 

Another important contribution to the efficiency comes from the dE/dx threshold cuts 
and the mass cuts. The inefficiency due to the dE/dx requirements were estimated sep- 
arately for positive and negative tracks using A — > pyr (+c.c.) decays. In this approach, 
protons were identified from the A peak and then the proton dE/dx selection efficiency 
was reconstructed as the ratio of the events without and with the dE/dx requirement. 
These efficiencies were determined as a function of p/M. The efficiency for each pt/M 
bin was corrected by reweighting the p/M distributions using Ariadne. The average 
efficiency of the dE/dx cuts for d{d) is 0.7 for px/M < 0.5. For larger momenta, the effi- 
ciency decreases due to the dE/dx > 2.5 mips cut. The signal extraction is not possible 
for pt/M > 0.7 due to a very small efficiency. For the low-momentum region px/M < 0.5, 
the efficiencies for negative tracks tend to be larger than for positive tracks. The dE/dx 
efficiency for p{p) is higher by 15% than that for d{d). 

Alternatively, the overall tracking and the dE/dx efficiency was calculated using the 
Ariadne MC model; consistent results with the approach discussed above were found. 

8 Systematic uncertainties 

The systematic uncertainties were evaluated by changing the selection and the analysis 
procedure. Only the largest contribution of each cut variation for the final invariant cross 
section is given below. The following sources of systematic uncertainties were studied: 

• efficiency of the track reconstruction and selection. The systematic uncertainty on the 
tracking efficiency for p, p, d was ±2%. This systematic uncertainty was found after 
variations of the track-quality cuts. For d, the systematic uncertainty, ±5%, includes 
both the effect of track-quality-cut variations and the reduction in e{d) when the linear 
model for the d absorption was used (see Section mi) : 

• efficiency due to the dE/dx selection. This systematic uncertainty was estimated by 
varying the cut dE/dx > 2.5 mips within the dE/dx resolution and by using the MC 
simulation. This systematic uncertainty was ±5%. For the lowest px/M bin, the 
uncertainty was ±10%; 
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• variations in the particle yields associated with the signal extraction: 

— the number of d{d) were reconstructed using a Gaussian fit to the DC A distribu- 
tion with a first-order polynomial for the background description; 

— the region used to determine the background for the side-band background sub- 
traction was reduced to 1.5 <| DC A \ < 3.5 cm; 

— the DC A cut for the side-band background subtraction was varied within its res- 
olution of ±0.1 cm; 

— for the side-band background subtraction, the background shape was taken from 
the MC (without d at the generator level); 

— the cut on AZ was varied by ±0.2 cm; 

These variations lowered the production yields by 5.0% for p, 2.2% for p, 26.0% for d 
and 6.1% for d. The largest effect originates from the conservative treatment of the 
shape of the DC A background. The upper systematic error was below 1% for p, p 
and d, and 11% for d. 

• the background contribution under the Z^tx peak for d events was assumed to be due 
to beam-gas interactions and, therefore, it was subtracted (—4% contribution for p, p, 
d and —17% contribution for d); 

• the correction for A decays applied for the p{p) sample was changed by ±10% (see 
Section [Q.ip . The size of this uncertainty, which is similar to that in other publica- 
tions [4,9], was determined by the uncertainty on the strangeness suppression factor 
in the Ariadne model; 

• variations of the DIS-selection cuts. The cut on the energy of the scattered electron 
was increased to 10 GeV, and the lower cut on the 6 distribution was tightened to 
40 GeV. The cut on Z^tx was varied by ±5 cm. The cut on the number of primary 
tracks was increased from three to four. These variations led to changes of for 
P, -I'a^/o for p, for d and ^.1^,3% for d. Variations of the cuts on ye and yjB 
distributions showed a negligible effect. 

The overall systematic uncertainty was determined by adding the above uncertainties in 
quadrature. The largest experimental uncertainty was due to the uncertainties on the 
tracking efficiency and the signal extraction. 

9 Results 

9.1 Production cross sections and B2 

For each particle type i, the invariant differential cross section can be calculated from the 
rapidity range Ay and the transverse momentum p^^i of a corresponding particle through 

j^^^a^ ^ J. 1 A^^ 

atot d(p./M,)3 A^Dis 27r{pT,i/Mi)Ay A{pT,i/Mi) ' 
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where the subscript i denotes a p{p) or a d{d), Ni is the particle yield in each pT,i/Mi 
bin after the correction for the tracking efficiencies and the particle selection and iVDis = 
2.59 X 10^ is the number of DIS events used in the analysis. For the present measurement, 
Ay = 0.8 and A{pT,i/Mi) = 0.1 are the bin sizes. For comparisons with other experiments, 
the p{p) rate was corrected for the decay products of A. A correction factor of 0.79 was 
estimated from the Ariadne simulation which gives an adequate description of Kg and 
A production [38]. 

The invariant differential cross sections as a function oi pt/M for p{p) and d{d) are shown 
in Fig. [7] and given in Tables [1] and [2l The d{d) invariant cross section is smaller by 
approximately three orders of magnitude than that of p{p). These cross sections were 
used to extract the coalescence parameter B2 as discussed in Section |2j The parameter 
B2 is shown in Fig. [8] and listed in Tables [3] and HI For d, B2 tends to be higher than for 
d, especially at low px/M. The value of B2 for d is in agreement with the measurements 
in photoproduction [9], but larger than that observed in e"'"e~ annihilation at the Z 
resonance [4]. The measured B2 is also significantly larger than that observed in heavy- 
ion collisions [11]. 

The events containing at least one p{p) or d{d) were analysed in the Breit frame [39]. The 
number of events with p{p) in the current region of the Breit frame was about 2.5% of 
the total number of observed events with p{p). In this region, neither d nor J was found. 
Since the current region of the Breit frame is analogous to a single hemisphere of e~^e~ , 
the observation of d{d) reported in this paper is not in contradiction with the low d rate 
observed in e~^e~ [2-4]. 

9.2 Production ratios 

The detector-corrected d/p and d/p ratios as a function oi Pt/M are shown in Fig. [9](a) 
and listed in Tables [3] and |H For the antiparticle ratio, there is a good agreement with 
the HI published data for photoproduction [9], as well as with pp data [8]. A similar d/p 
ratio was also observed in hadronic T(IS') and T(2S') decays [2]. 

The d/d and p/p ratios as a function oipx/M are shown in Fig. [9]^b) and listed in Table O 
The p/p ratio is consistent with unity, as expected from hadronisation of quark and gluon 
jets. The dominant uncertainty on the ratio is due to systematic effects associated with 
the track selection and reconstruction. 

The production rate of d is higher than that of d, especially at low pt- Under the 
assumption that secondary interactions do not produce an enhancement at DC A = for 
the d case, the result would indicate that the relation between d/d and {p/p)"^ expected 
from the coalescence model does not hold in the central fragmentation region of ep DIS 
collisions. 

For collisions involving incoming baryon beams, there are several models [40, 41] that 
predict baryon-antibaryon production asymmetry in the central rapidity region. A p — p 
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asymmetry in proton- induced reactions is predicted to be as high as 7% [41]. Given the 
experimental uncertainty, this measurement is not sensitive to the expected small p — p 
asymmetry. 

In heavy- ion collisions, the dtod production ratio is expected to be smaller than unity [42]. 
A recent measurement at RHIC [11] indicated a lower production rate of d compared to 
that of d. The average value of the ratio d/d = 0.47 ± 0.03 was compatible with the 
square of the p/p = 0.73 ± 0.01 ratio. Assuming the same size of the production volume 
for baryons and antibaryons, this RHIC result is consistent with the coalescence model. A 
similar conclusion was obtained earlier in fixed-target pp [8] and pA [7] experiments. For 
e~^e~ collisions, the d yield is compatible with that of d within the large uncertainties [4,5]. 

10 Summary 

The first observation of d{d) in ep collisions in the DIS regime at HERA is presented. The 
production rate of d{d) is smaller than that for p{p) by three orders of magnitude, which 
is in broad agreement with other experiments. 

The production of d{d) was studied in terms of the coalescence model. The coalescence 
parameter is in agreement with the measurements in photoproduction at HERA. However, 
it is larger than that measured in e^e~ annihilation at the Z resonance. 

The production rate oip is consistent with that of p in the kinematic range 0.3 < px/M < 
0.7. Due to significant uncertainties, it is not possible to test models that predict a small 
baryon-antibaryon asymmetry in the central fragmentation region. 

For the same kinematic region, the production rate of d is higher than that for d. If 
the observed d are solely attributed to deuterons produced in primary ep collisions, the 
results would indicate that the coalescence model with the same source volume for d and 
d cannot fully explain the production of d{d) in DIS. 
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m„ 1 A/f 

PT 1 Ivi 




[Id/ <^tot)(i CTd/ (i[Pd/ -^y^d) l^XiU ) 


0.3 - 0.4 


1 o o 1 r\ r\ ^ — r-( 1 1 4 

1.33 ± 0.01_o;2i 


3.29 ± 0.43_i;24 


0.4 - 0.5 


1 o yi 1 r\ r\ 1 -l-0,16 

1.34 ± 0.01_o;];g 


1.37 ± 0.26_o;5{ 


0.5 - 0.6 


0.88 ± 0.01^!^;^^ 


1.16±0.28l[5i^ 


0.6 - 0.7 


0.38 ±0.01+°;°^ 





Table 1: The measured invariant cross sections for the production of p and d in 
DIS as a function of px/M. The statistical and systematic uncertainties are also 
listed. 



Pt/M 


(7p-/o^tot)dVp/d(pp/Mp-)3(xlO-2) 


(7j/atoi)dVj/d(pj/M,-)3(xlO-5) 


0.3 - 0.4 


1.59±0.0llJ^i^ 


0.77 ±0.15^°:?^ 


0.4-0.5 


1.21 ±0.01^°:°^ 


0.45 ± o.iii[5:[!? 


0.5 - 0.6 


0.86 ± o.oiij^:^!^ 


0.60 ± o.i9l[5:[!^ 


0.6 - 0.7 


0.35 ± o.oilj^:^!^ 





Table 2: The measured invariant cross sections for the production of p and d in 
DIS as a function of px/M . The statistical and systematic uncertainties are also 
listed. 



Pt/M 


R{d/p){xlO-^) 


52(rf)(10-2GeV2) 


0.3 - 0.4 


2.48 ± 0.33^?:^^ 


4.11 ±0.541};^^ 


0.4 - 0.5 


1.02 ± 0.19lJ^i^ 


1.68 ± 0.321°:^° 


0.5 - 0.6 


1.32 ±0.32+°!^ 


3.31 ± 0.801^5 


0.6 - 0.7 






0.3 - 0.7 


1.88 ± 0.20lJ^;^^ 


3.32 ± 0.34+};^^ 



Table 3: The measured d-to-p production ratio and the parameter B2 for d as 
a function of pt/M . The last row of the table shows the data in the full measured 
phase space. The statistical and systematic uncertainties are also listed. 



Pt/M 


i?(rf/p)(xl0-3) 


52(c/)(10-2GeV2) 


0.3 - 0.4 


0.48 ± 0.09lJ]:?^ 


0.67±0.13lJ]:|^ 


0.4 - 0.5 


0.37 ± o.o9lJ^:[!^ 


0.67±0.17lJ^:i^ 


0.5 - 0.6 


0.70 ± 0.22+°:°^ 


1.80±0.57l°J^ 


0.6 - 0.7 






0.3 - 0.7 


0.49 ± 0.07+°:°^ 


0.89 ±0.14+°:^^ 



Table 4: The measured d-to-p production ratio and the parameter B2 for d as 
a function of pt/M . The last row of the table shows the data in the full measured 
phase space. The statistical and systematic uncertainties are also listed. 
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Pt/M 


R{p/p) 


R{d/d) 


0.3 - 0.4 


1.19 ±0.01^!]:?^ 


0.23 ± 0.05t!]:|]^ 


0.4 - 0.5 


0.90 ±0.01+°:^° 


0.33 ±0.10+°:^? 


0.5 - 0.6 


0.97±0.0i;[JiJ 


0.52±0.2li°:i^ 


0.6-0.7 


0.92 ± 0.03^!^:^^ 




0.3-0.7 


1.05±0.0ltU;}^ 


0.31 ± 0.05+!^:^;^ 



Table 5: The measured p-to-p and d-to-d production ratios as a function of 
Pt/M . The last row of the table shows the data in the full measured phase space. 
The statistical and systematic uncertainties are also listed. 
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Figure 1: The dE/dx distributions as a function of the track momentum for (a) 
positive and (h) negative tracks. The DIS events were accepted by requiring at least 
one track with dE/dx > 2.5 mips (denoted by the dashed lines), |AZ| < 1 cm and 
\DCA\ < 0.5 cm. The lines show the most-probable energy loss calculated using the 
Bethe-Bloch formula for different particle species. 
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Figure 2: The mass spectra for (a) positive and (b) negative particles. Tracks 
are selected as for Figure [21 The mass distribution was calculated from the track 
momenta and the dE/dx. The arrows indicate the cuts applied for the selection of 
candidates. 



19 



0) 

■Z 40000 

c 

LU 




AZ (cm) 



2 4 6 

AZ (cm) 



Figure 3: T/ie distributions of AZ , the distance of the Z -component of the track 
helix to Zvtx for: (a)-(h) particles and (c)-(d) antiparticles, as indicated in the 
figure. The p, p, d and d candidates were identified using the dE/dx mass cuts 
(see text). The arrows indicate the applied cuts. 
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Figure 4: The distributions of the distance of closest approach, DCA, for: (a)- 
(h) particles and (c)-(d) antiparticles. The DCA are shown after the cut \ AZ |< 
2(1) cm as discussed in the text. The arrows indicate the signal region for the side- 
hand background subtraction. The dashed lines show the fitted background level. 
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Figure 5: The distributions of the number of events with at least one d(d) andp(p) 
candidate normalised to unity as a function of: (a)-(d) DIS kinematic variables and 
(e) rapidity y. The points for d and d are slightly shifted horizontally for clarity. 
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Figure 6: The Z^t^ distributions for: (a)-(h) particles and (c)-(d) antiparticles, 
as indicated in the figure. The solid lines show the fit using a Gaussian distribution 
with a first-order polynomial function for the background description. The dashed 
line shows the fitted background. The arrows indicate the cuts applied for the final 
selection. 
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Figure 7: The invariant differential cross sections for p{p) and d{d) produced in 
DIS ep collisions as a function of px/M. The inner error bars show the statistical 
uncertainties, the outer ones show statistical and systematic uncertainties added in 
quadrature. For clarity, the points for particles and antiparticles are slightly shifted 
horizontally with respect to the corresponding pt / M . 
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Figure 8: The px/M dependence of the parameter B2 for d and d produced in DIS 
ep collisions and in photoproduction [9]. The inner error bars show the statistical 
uncertainties, the outer ones show statistical and systematic uncertainties added in 
quadrature. For clarity, the points for particles and antiparticles are slightly shifted 
horizontally with respect to the corresponding pr / M . 
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Figure 9: (a) d/p and d/p production ratios as a function of px/M compared 
to the HI photoproduction results [9]. (h) the d/d andp/p production ratios as a 
function ofpr/M. The inner error bars show the statistical uncertainties, the outer 
ones show statistical and systematic uncertainties added in quadrature. The points 
in (a) are slightly shifted horizontally for clarity. 
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